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Abstract. We present measurements of the mid-infrared (MIR) emission from SN 1987A, made using the Infrared 
Space Observatory (ISO) 11 years after outburst. They are the only late epoch detections of this source in the 
thermal IR regime. The position of the source, determined from an offset to an IR-emitting star, suggests that the 
emission is associated with SN 1987A or its extended supernova remnant (SNR). A predominantly circumstellar 
origin is however suggested by the size and orientation of the IR-emitting region, which is comparable with the 
extension of the inner ring seen with the Hubble Space Telescope (HST). The emission is most probably from 
coUisionally-heated circumstellar grains embedded in shocked gas downstream of the blast wave. The MIR extent 
is consistent with the hypothesis that the blast wave was propagating into material of moderate density interior 
to the thick inner ring at the epoch of the ISOCAM observations. 
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1. Introduction 

Supernovae are thought to play an important role in the 
life cycle of interstellar dust grains. On the one hand, 
shocks driven by supernovae into the interstellar medium 
(ISM) are considered to be responsible for the destruc- 
tion of interstellar grains (see e.g. ShuU, 119771 Draine & 
Salpeter. [T979bl Tielens et aL lTOMl Jones et aL lTOMjl . 
On the other hand, supernovae have also been proposed as 
sources of interstellar dust (e.g. Hoyle & Wickramasinghe, 
[TOTDl Dwek & Scalo. lTMTji . 

SN 1987A in the Large Magellanic Cloud (LMC) was 
the first supernova in which the condensation of dust in 
the central metal-rich part of the expanding ejecta has 
been observed. The appearance of a continuum emission 
in the IR in the spectrum of SN 1987A 614 days after 
outburst is attributed to grains formed after at least 500 
days (e.g. Wooden. 11997(1 . The infrared emission from the 
condensates was (maybe except at the onset of dust for- 
mation; Roche et al., I1993|l optically thick at the early 
epochs when the IR-emission was still bright enough to 
be observable from ground. Therefore, it was not possi- 
ble to get clear information about their properties and in 
particular their mass. 
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Observations with the HST have revealed a bipolar 
ring structure consisting of a thick inner ring centred at 
the position of the supernova, and two thin outer rings of 
larger radius at each side (Burrows, T995). Radio and X- 
ray observations suggest that, after an initial period of free 
expansion, the blast wave started to propagate through a 
region of moderate density with njj ~ 100 cni^^ interior 



to the thick inner ring (Chevalier & Dwarkadas, I1995|l 
Chevalier & Dwarkadas suggested this region to be an 
HII region, illuminated by the blue supergiant supernova 
progenitor and composed of wind material ejected during 
the prior red supergiant phase. 

If the HII region contains dust, extended thermal emis- 
sion from coUisionally heated grains embedded in the 
shocked gas downstream of the blast wave is to be ex- 
pected. The presence of circumstellar grains was already 
invoked at early epochs to account for observations of 
a prompt echo effect in the mid-infrared (Rank et al., 
119881 Chabalev et al.. I1989|l) . It has also been suggested 
that circumstellar grains can account for observations af- 
ter 580 days showing an extension of 1'.'5 at 9 /im (Roche 
et al.. ll998|l . Light echos in the optical have further been 
used to reconstruct a three dimensional model of the larger 
scale structures around SN 1987A (Crotts et al.. 11995(1 . 

The ISO mission (Kessler et al., 19965 gave a unique 
opportunity to probe the early evolution of SN 1987A in 
the IR prior to the envelopment of the thick inner ring 
by the blast wave. In this paper we present sensitive mea- 
surements made with the ISOCAM instrument (Cesarsky 
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et al.. I1996|l on board ISO 11 years after outburst. Thus 
far, these measurements constitute the only late epoch de- 
tection of the source in the IR regime. Preliminary results 
were presented by Tuffs et al. p998|l . Here we demon- 
strate that the detected IR emission is most probably of 
circumstellar origin. 

The paper is structured as follows: In Sect. 2 we 
present the ISOCAM measurements and describe the re- 
duction and calibration of the data. In Sect. 3 we derive 
the flux density, position, orientation and size of the mid- 
infrared (MIR) source, showing that the IR emission from 
the remnant is resolved, with an extent comparable to that 
of the thick inner ring. These results are discussed in Sect. 
4. A summary of the paper is given in Sect. 5. 

Throughout the paper we adopt a distance of 51 kpc 
to the supernova as was used to analyse the X-ray ob- 
servations by Hasinger et al. H1995|l . This is close to the 
distance of 51.2 ±3.1 kpc derived by Panagia et al. fT991) 
but slightly larger than the values given in later publica- 
tions (Gould, Cnni Sonneborn et aL. ITWl Gould. 11998 ). 

2. Observations 

The observations were made with the ISOCAM instru- 
ment in the LW2, LWIO and LW3 broad-band filters, 
whose pass bands are centred on 6.75, 12 and 14.3 /im re- 
spectively. An observation log summarizing the principal 
observational parameters is given in table □ 

The measurements in LWIO and LW3 were optimized 
for obtaining sensitive photometry and structural infor- 
mation of a faint compact target. A pixel field of view 
(PFOV) of 3" was chosen as for the anticipated back- 
grounds it gave the best trade off between angular res- 
olution and detector illumination. We used a 3 x 3 raster 
with raster interval 10" x 10" to sample the sky at inter- 
vals of 1". This comfortably oversampled the point spread 
function (psf), which has a FWHM of ~ 4" at 10 fim for 
the 3" PFOV. 

The raster interval of 10" was the minimum interval 
that provides a sky sampling for which each view of the 
target in a given pixel is immediately preceded or followed 
by a view of the underlying background. This optimized 
the knowledge of the relative response of the detector be- 
tween source and background, which in turn could be used 
to minimize the effect of residual uncertainties in the flat- 
field response of the 32 x 32 pixel detector on the de- 
termination of source structure. The observations were 
separated in epoch over the period September 1997 to 
February 1998. The measurement in the 12 fim band was 
a follow up observation to the detection in the 14.3 /im, 
while the measurement in the 6.75 /im band was part of an 
independent program of large-scale mapping of the LMC 
led by L. Vigroux. 

2.1. Calibration and Data reduction 

The data were first corrected for dark current using 
the model of the Cam Interactive Analysis (CIA) pack- 



Table 1. Observation log 



Filter 


LW2 


LWIO 


LW3 


TDT 


667017 


811022 


750019 


observer 


L. Vigroux 


R.J. Tuffs 


R.J. Tuffs 


Arof [/xm] 


6.75 


12. 


14.3 


filter range [fim] 


5.00-8.50 


8.00-15.00 


12.0-18. 


observation date 


13.9.97 


5.12.97 


3.2.98 


day after outburst^ 


3855.2 


3998.7 


3937.7 


raster 


2x2 


3x3 


3x3 


pixel size 


6" 


3" 


3" 


step size 


24" 


10" 


10" 


read out interval tint 


2.10 s 


5.04 s 


5.04 s 


reads per pointing ni 


31 


19 


19 



" The time of the outburst is taken to be February 23.316, 
1987 (Arnett et aL. nM^ . 



age (Ott et aL. ITOMI version April 2000). This was fol- 
lowed by glitch removal using the Multi resolution Median 
Transform^ also provided by CIA. The data were then in- 
spected and residual events and strong longer-lived glitch- 
induced distortions were removed manually. To allow an 
accurate determination of the illumination seen by each 
single pixel during the raster, the data were then corrected 
for the transient response behavior of the detector. To this 
end we developed and implemented a routine based on 
the drift model of the ISOPHOT-S detector on board ISO 
f Schubert. I1995|l . The routine derived the time constants 
and jump factors determining the drift behavior of the 
individual detector pixels from the measured data. 

To create a map showing the background structure 
underlying the sources we applied a flat-field correction, 
which was optimized for the smooth variation of the back- 
ground. The flat-field was derived from a map, which 
had already been corrected using a library flat provided 
by CIA and which had been spatial filtered to reduce 
the noise. The fiat-field-corrected data were then binned 
(without interpolation) on a grid with a resolution of 1". 
For visualization the data were further filtered by setting 
pixel values differing more than one sigma from their sur- 
roundings to the average value of the surrounding pixels 
(the filter was not applied for obvious discrete sources). 
Then map pixels without any data were set to the av- 
erage values of the neighbouring pixels. Finally, we con- 
volved the image with a Gauss function with a variance 
of 0'.'6. The resulting maps are shown in Fig.^ 

The quantitative analysis of the individual compact 
sources seen in the ISOCAM maps was done directly from 
the observed brightnesses as a function of sky position, 
without use of maps, as described in Sect.|21 This required 
a different calibration procedure compared to that for the 
maps. For the observations at 12 and 14.3 //m we made 
use of the almost flat emission around the two sources 
and normalized all pixel values at individual raster po- 
sitions with the corresponding values of the background. 
This was done by making a linear fit to the variation of 
(responsivity drift-corrected) illumination in each given 
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a (2000) 5h a (2000) 5h 



Fig. 1. Region around the SN 1987A in the MIR, as seen at 12 {LWIQ) and 14.3 /xm {LWi). The brightness 
increases with darkness and ranges from 0.215 to 0.306 mJy/["]^ {LWIQ) or from 0.356 to 0.426 mJy/["]^ {LWi) 
respectively. The source close to the middle is emission from the remnant of SN 1987A. The IR-emission at the border 
of the image is the reference star, that we used to derive the absolute coordinates of the IR-emission from the remnant 
of SN 1987A. The bright emission in the north is an arm of 30 Doradus. The two individual sources show different 
colours, suggesting a lower temperature for the IR-emission from SN 1987A than for the star. 



detector pixel with spacecraft raster position, leaving out 
raster positions where the pixel was viewing the source. 
In the case of the measurement at 6.75 /im, we used the 
library flat provided in CIA, as the raster pattern was 
not suitable for this method. For the measurement at 12 
and at 14.3 /im the sky coordinates of each detector pixel 
at each raster pointing were evaluated from the actual 
(post posteri) pointing position of the satellite using the 
median of all the data used. To derive the pixel coordi- 
nates of the observation at 6.75 /im we made use of the 
astronomical information provided by the CIA-program. 
Corrections were applied for the astrometrical distortions 
in the ISOCAM field (Aussel, induced by the field 

lens using the polynomial correction coefficients for the 
different filters as given in CIA. 

2.2. The ISOCAM maps 

The final images show two clear individual sources. The 
source close to the middle corresponds to emission from 
the remnant of SN 1987A, whereas the source near the 
edge is an IR-emitting star. The two sources show differ- 
ent colours, indicating a higher temperature for the star. 
The serendipitous detection of the star provided a unique 
possibility to determine an accurate position for the cen- 
tral MIR source near SN 1987A (see Sect.ESl- The bright 
ridge of emission to the NNW of SN 1987A is an arm of 
the 30 Doradus nebula. Although backgrounds are dom- 
inated by smooth emission from zodiacal light, there is 
some evidence for faint cirrus structure at the 1% level in 



the southern part of the field. At both wavelengths the IR 
emission near the map centre is offset several arcsec from 
the nominal position of SN 1987A. These offsets are partly 
due to the inaccuracy of the position of the satellite but 
are mainly attributable to a displacement of the field lens 
of the camera (Blommaert & Cesarskv. ll998|l . 

3. Parameters of the discrete sources 

Various parameters of the discrete sources, expressed as 
the vector a, were derived using a non-linear x^-fii, where 
we compared the measured data f{xi) and uncertainties 
Gi at the sky positions Xi around the individual sources 
with a model image B{xi, a): 

2. N {B{x„a) - f{x,)f 
X{a)^2^ ^ • (1) 

The model image B{x,a) is in general the convolution 
of the theoretical psf corresponding to the telescope op- 
tics and the spectral response of the used filter band pass 
with the source function (taken to be a Gaussian), on a 
certain background. For a symmetrical Gaussian fit a is 
[Fq, Xo,aQ, Ijj] and the model image given by: 

B{x,a) = + J Jp{x-x') 

"2^^"P(^ ' 

where Jf{x) is the theoretical broad-band psf (see ap- 
pendix^, Fq the flux density, Xq the coordinates, ctq the 
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position In Y 




position in Y 



Fig. 2. Best fits of the theoretical image (solid line) to the measured brightnesses (diamonds) of the remnant of 
SN 1987A at 12 /xm (LWIO) and 14.3 /xm {LW3). They are shown as scans along lines (at fixed Z in the coordinate 
system of the satellite) through a subimage centred on SN 1987A. The dotted vertical lines mark the beginning and 
the end of the scans. To calculate the psf for each filter the spectral shape was assumed to be Fx oc A"*^ with a 
spectral index k = 1.08. The FWHM of the circular Gaussian source is 1'.'50 (LWIO) and 1'.'13 {LW3). 



variance of the Gaussian source and lu the background. 
The expression n was minimized to derive Fq, xq, Ijj (at 
6.75 /im only) and aq (only at 12 and 14.3 /xm). 

To achieve a higher accuracy, in particular for the de- 
termination of the size of the emitting region, we took into 
account the dependence of beam profile on source colour. 
This was done by averaging simulated beam shapes ob- 
tained at a series of wavelengths across the broad filter 
pass bands, weighted according to a power law source 
spectrum -FaA*^ = const, with a spectral index k (see 
appendix^). For the star the spectral index was always 
taken to be fc = 4 appropriate to a Rayleigh- Jeans spec- 
trum. For the source associated with the supernova k was 
estimated from the finally determined colour corrected 
(see e.g. Blommaert, 1998) flux densities to be 1.08±0.32. 
The model image was calculated on a grid with a resolu- 
tion of 1/3". In the fitting procedure the predicted fluxes 
at the observed positions were found using linear interpo- 
lation. 

The true uncertainties ai in each individual measure- 
ment for a detector pixel at a given spacecraft raster point- 
ing direction arc difficult to estimate and may be partly 
affected by systematic uncertainties. This particularly ap- 
plies to measurements viewing the discrete sources, where 
the illumination and detector response function may vary 
with position over the field of view of the individual detec- 
tor pixels. We therefore assumed in nearly all model cal- 
culations that the uncertainties of all measurements are 
equal. This is justified by the fact that the distribution 
of the residuals was found to be normal. The only excep- 
tion to this was made for the measurement of the star at 
6.75 /im, where the uncertainties of the high source fluxes 
were mainly attributed to a particularly strong transient 
behavior. 



The absolute uncertainties <Ti of the measured data 
were chosen such that the value of the minimum (Eq.[T]l 
corresponds to a confidence level of 90%. In the case of the 
measurements at 12 and at 14.3 /im we made a special 
fit letting Fq, Xq, cq and k as free variables. The joint 
uncertainties of v parameters were estimated by finding 



the range within which the function ^xt — xt 



did 



not exceed the value corresponding to a confidence level 
of one sigma, letting all other parameters as free variables. 
Thus, apart from the jointly determined uncertainties in 
the coordinates, which were estimated from Ax| = 2.30, 
the uncertainties in individual parameters correspond to 
Ax? = 1 (see e.g. Press et aL. nM^ . 

3.1. Fits with a symmetrical Gaussian as IR-source 

In the case of the measurement at 6.75 /im the parameters 
varied were -Fq, the coordinates Xq and the background 
brightness /y. We used a model image corresponding to 
an ideal point source (crq = 0) assuming a spectrum with 
fc 1 for the IR emission from the supernova remnant. 

At 12 and at 14.3 /im the flux density and the coor- 
dinates were optimized together with the extension of the 
source. The variables varied were therefore Fq, Xq and 
(Tq . Because of the normalization of the pixel data to the 
background we held /u fixed. The uncertainties in Fq and 
Xq were estimated for a fixed (most probable) source size. 
The best fits to the measured data with LWIO and LW3 
are shown in Fig. |5| 

The surface brightness of the background at these two 
wavelengths was estimated from the images calculated as 
described in Sect. 2.1, but without filtering and convolu- 
tion. To derive the surface brightness Itj at the position 
of the supernova from the maps we took the average value 
within two circular areas of radius of 5" offset 15" to the 



Fischera et al.: Infrared emission towards SN 1987A 11 years after outburst 5 



Table 2. Spherical Gauss fits to the IR emission associated with SN 1987A and the reference star 





SN 1987A 


reference star 


Filter 

ADU7(mJy/pixel)'' 
beam size [FWHM ["])" 


LW2 
6.75 

2.32 ± 3.3% 
6.887 ± 0.002 


LWIO 
12. 

4.23 ± 3.9% 
4.195 ± 0.025 


LW3 
14.3 

1.96 ± 4.8% 
5.109 ±0.009 


LW2 
6.75 

2.32 ± 3.3% 
6.222 ± 0.001 


LWIO 
12. 

4.23 ± 3.9% 
3.945 ± 0.026 


LW3 
14.3 

1.96 ± 4.8% 
4.955 ± 0.008 


spectral index fc'' 

derived source size {FWHM ["])" 
xl (^) 

Uncertainty of pixel values a/Iv 
Fq/Iv [['f] 

background 7u [mJy/["]^] 
source flux density Fq [mjy]* 
source flux density Fq"''" [mjy]'' 


1. 

0.76 (51) 
3.6% 

6.54 ± 25% 
0.0759 ± 3.3% 
0.50 ± 25% 
0.43 ± 25% 


1.08 

l.dU_(j 24 

0.82 (94) 
1.9% 

10.62 ± 1.8% 
0.226 ± 4.0% 
2.40 ± 4.3% 
2.64 ± 4.3% 


1.08 

< 2.3 (3 cr) 
0.83 (127) 
1.7% 

7.608 ± 2.6% 
0.370 ± 4.8% 
2.82 ± 5.5% 
2.90 ± 5.5% 


4 

0.77 (60) 
12.%^ 

116.2 ±5.5% 
0.0717 ± 3.6% 
8.33 ± 6.2% 
8.43 ± 6.2% 


4 

0.69 (41) 
2.5% 

16.11 ±3.1% 
0.228 ±6.1% 
3.67 ± 5.0% 
3.73 ± 5.0% 


4 

< 2.4(3cr) 
0.78 (75) 
1.5% 

4.879 ± 5.0% 
0.365 ± 5.0% 
1.78 ± 7.0% 
1.69 ± 7.0% 



" Analog to Digital Units 
Siebenmorgen et al.. 119991 

Average FWHM of cuts in Z and Y through the maximum of the theoretical psf. 
Uncertainties correspond to the difference of these two values. 

Spectral index assumed in calculating the theoretical psf (see text). 
° FWHM of the symmetrical Gaussian model assumed for the source. Uncertainties are given as 1 a. 
^ Median of the uncertainties. 
^ Flux density without colour correction. 

Colour-corrected according to a spectrum Fx oc X^^Bx{T) with T = 198.7 K (supernova) and T = 564.1 K (star). 



east and to the west of SN 1987A. At both 12 and 14.3 
the difference between the two values is less than 1%. For 
the background brightnesses at the position of the star 
we took the average of a number of single measurements 
around the star; here the estimated uncertainties in the 
background are 1.3% (LWIO) and 0.9% {LWS). 

The derived values for the source flux and size, and 
the background brightnesses, are summarized in table 12 
together with the assumed value for the spectral index k 
used in the determination of source size. The table also 
gives the source flux densities normalized to the surface 
brightness of the background. The uncertainties quoted 
for the normalized flux densities at 12 and 14.3 /im are 
statistical uncertainties derived from the x^-fit- 

The IR spectrum of the supernova remnant can be 
approximated by a simple black body spectrum with 
T = 288 ± 16 K (x^i„ = 0.84) or a modified black body 
spectrum Fa oc X~'^Bx{T) with T = 198.7i?J K where 
a good fit is achieved (Xmin = 0.02). Assuming the same 
modified spectrum for the star we get a temperature of 
T = 564.ltlll K. The fit is again very good (x^^in = 
1.2 X 10^'^). A simple black body for the star seems to be 
unhkely (xii„ w 5.3). 

3.1.1. Sizes of the IR-sources 

The probability distributions for the extension of the IR- 
emitting regions associated with SN 1987A were derived 
at 12 and at 14.3 /im from a x^"fit using a symmetrical 
Gaussian as source function (see appendix . These are 
shown in Fig. |3| The different probability distributions 
correspond to different colour indices used in the calcula- 
tion of the theoretical psf. The distributions plotted with 



thick solid lines have been derived for the most probable 
colour of the source, a power law with k ~ 1.08. As can 
be seen in the figure, the probability distributions are only 
slightly affected by the uncertainty in the spectral index k 
for the MIR emission associated with the supernova. The 
probability distributions corresponding to the upper and 
lower limits for k (0.56 and 1.59, respectively, at 90% con- 
fidence level) are plotted as thin solid lines. The dashed 
and dotted curves peaking near 1'.'9 are probability dis- 
tributions for the measurements at 12 /im using as psf 
the image of a calibration star and the theoretical psf for 
a spectral index k — A, respectively. Since the spectrum 
of the star should be a power law with = 4, the simi- 
larity of these two curves demonstrates that the effect of 
interchanging the theoretical psf with the measured psf 
is small. The larger apparent source sizes calculated from 
the psfs with fc = 4 arise from the decrease in FWHM 
of the psf with increasing k compared with the true value 
k = 1.08 ± 0.32 for the MIR emission associated with 
SN 1987A. 

For the LWIQ observation we obtained a FWHM of 
1.50^0 24 arcsec for the extension of the source. This size 
is in agreement with the extensions of the projected in- 
ner ring, seen in the HST picture. The distribution we 
derived from the LW3 observation is broader, and also 
has a higher probability for small extensions. This might 
be due to the strong glitches causing long time distor- 
tions of the detector response of the observation (see also 
Sect. l3.2|l . From this observation we were only able to ob- 
tain an upper limit for the source size of 2'.' 3 at 3a level. 
Table 121 gives the maximum and minimum sizes of the IR- 
emitting region for certain confidence levels derived from 
the probability distributions for k = 1.08. As shown in 
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appendix IB . 21 the theoretical psf for the LWIO filter is a 
sufficiently accurate representation of the real psf to re- 
solve sources with a FWHM greater than ~ 0'.'6. If we 
take the measurement at 12 /im the source is extended 
with a certainty of more than 99%. 
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Fig. 3. Derived probability distributions of the size s of 
the theoretical image (see Eq. EH) as a function of the 
FWHM of the IR source for the two filter measurements 
LWIO (12 fim) and LW3 (14.3 iim) of SN 1987A. The 
emitting region is assumed to be a circular Gaussian. 
Also shown is the effect of different source spectra of form 
Fx (X X^'^ with spectral index k. The probability distribu- 
tions shown as thick lines correspond to the most probable 
spectral index k = 1.08 as derived from the flux densities 
at 12 and 14.3 ^m. The distributions corresponding to the 
upper and lower limits for k (0.56 and 1.59, respectively, 
at 90% confidence level) are plotted as thin solid lines. 
The probability distributions shown as dashed and dotted 
lines are size measurements of the IR emission associated 
with the supernova at 12 /xm using as psfs an image of 
a calibration star (HIC78527) and a theoretical psf with 
k — A, respectively. 



Table 3. Sizes of the IR source associated with SN 1987A 



conf. level 


68 % 90 % 95 % 99 % 99.9 % 




Maximum size of the source {FW H M) 


LW3 
LWIO 


1'.'46 1'.'79 1'.'92 2'.' 15 2'.' 40 
I'.'eO 1'.'76 1'.'83 1'.'95 2'.' 08 




Minimum size of the source {FW H M) 


LW3 
LWIO 


I'.'OS 0'.'62 0'.'45 0'.'21 0'.'07 
1'.'40 1'.'19 1'.'09 0'.'87 0'.'56 



Some evidence was found for extended emission asso- 
ciated with the star at 12 /im (Table ISJ. This may in part 
be due to the strong edge position and the course sam- 
pling of the star. But it cannot be excluded from the IR 
observations that the star itself is extended. 



3.2. Fits with an elliptical Gauss function as IR-source 

In addition to the extension we also examined the orien- 
tation angle of the IR-source. Herefore we described the 
source through a two dimensional elliptical Gauss function 
with variances crj and ay and a rotation angle which 
was taken to be the angle between the negative Z-axis 
in the satellite coordinate system and the y-axis of the 
rotated Gaussian, measured counter clockwise, a is then 
[Fq, a;o, ctj, (Ty, t?, lu] and the model image given by: 



Fn 



27r crjCTy 



X exp 



2a| 



2a| 



(3) 



where x{d, x, Xq) — {x{-d, x, Xq), y{-d, x, Xq)) are the coor- 
dinates of the rotated Gaussian: 



x{'d,x,XQ) = 



cos('!9) sin(i9) 
- sin(z9) cos(i?) 



(X - Xq). 



(4) 



As before we used theoretical psfs for the two filters LWIO 
and LW3 corresponding to A; = 1.08. For each rotation an- 
gle 1} we derived best fit variances through the minimiza- 
tion of x^(o) (Eq.^, leaving Fq and Xq as free parameters 
and lu fixed. The position angle of the elliptical Gaussian 
is taken to be the angle of the y-axis to the northern di- 
rection, measured again counter clockwise. The variation 
of the minimum Ax^ with the position angle for the two 
filters is shown in Fig. 01 together with the corresponding 
extensions given as FWHM. The extensions and the posi- 
tion angle are compared with the orientation and the sizes 
of the minor and the major axis of the projected inner ring 
of SN 1987A, seen with the WFPC2 of the HST (Burrows 
et al., 1995). The source sizes of the best fit are shown 



in Fig. 5(b) overlaid on an WFPC2-image, where we have 
centred the ellipses representing the fitted Gaussian source 
model at the position of the SN 1987A. 

The variation in Ax^ is low in the case of the mea- 
surement with LWIO. But the coincidence of the derived 
parameters of the extensions and the rotation angle of the 
two dimensional Gauss function with the parameters of 
the inner ring are remarkable. The derived region of the 
IR source at 14.3 /im on the other hand seems to have 
an extension with the dimension of the inner ring only in 
one direction. In the other direction, the shown extension 
is the minimum we allowed to avoid inaccuracies in the 
modeling. We think that this strong asymmetry is artifi- 
cial and an effect due to strong glitches at the beginning 
of the LW3 observation. 

3.3. Position of the individual sources 

As an independent property of the IR emission from the 
remnant of SN 1987A we derived the absolute coordinates 
relative to the position of the IR-cmitting star visible in 
the two mosaics at 12 and even at 14.3 /im (Fig.^. The 
position of the IR-emitting star itself was not known so 
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Fig. 4. The top panel shows the variation of Ax^ with the position angle of a two dimensional Gaussian as source 
function of the IR emission at 12 /j,m (LWIO) and 14.3 fim (LWS). Below the FWHM of the two directions of the two 
dimensional Gaussian (solid and dashed line) are shown as function of the position angle. The extensions are compared 
with the lengths of the two axes of the elliptical projection (dotted lines) of the thick inner ring of SN 1987A, seen with 
the WFPC2 of the HST (Burrows et al.. I1995|l . The position angle of the major axis of the ring projection (81.2°) is 
shown in the top panels as the vertical dashed-dotted lines. The given confidence levels belong to the minimum Ax^. 



that we first had to determine the absolute coordinates of 
this star (see appendix Q. 

We derived the following offset positions of the IR 
emission from the position of the supernova, where we 
have taken the coordinates given by Reynolds et al. l|TM5|l : 



Aa 
AJ 

and 

Aa 
A,5 



-0'.'04±0'.'24stat. ±0'.'42 
-0'.'48 ± 0'.'24stat. ± 0'.'42 



sys., 
sys. 



-0'.'77±0'.'15.tat. ±0'.'42 
-0'.'09±0'.'15.tat. ±0'.'42 



sys., 
sys.- 



(LWS) 



(LWIO). 



(5) 



(6) 



The systematic uncertainty is not very well known. It is 
simply assumed that the systematic uncertainty in indi- 
vidual directions is given by the accuracy of the pixel po- 
sitions in the inner ISOCAM-field after distortion correc- 
tion, which should be accurate to within a tenth of a pixel 
(Aussel. I1998|l . The combined uncertainty of the coordi- 
nate pair is therefore taken to be ~ 0'.'42. Due to the strong 
edge position of the star in particular in the LW10-?ie\d 
and the fact, that at one raster point the dead line of 
the LW camera of the LWIO observation was lying at the 
sky position of the star, the systematic uncertainty of the 
offset position at 12 /im may be larger. 

The absolute coordinates of the IR-emission from the 
remnant of SN 1987A at 12 and 14.3 /ini with the cor- 
responding uncertainties are shown in Fig. 



5(a) 



rela- 
tion to the WFPC2-image. Fig. 5(b) shows a visualization 
of the fits to the size of the emitting region. As seen in 



Fig. 5(a) the IR emission measured with ISOCAM should 
be from of the remnant of SN 1987A. Taking the joint 



probability distribution of both measurements, the most 
probable location for the IR counterpart is to the north 
east of the supernova, which is also the position of the 
maxima in the X-ray and radio surface brightness distri- 
butions (Burrows et al. . 1200011 . But it is also possible that 
the centre of the IR emission lies at the position of the su- 
pernova so that the outer boundary of the emitting region 
could indeed coincide with the inner ring as seen in Fig. 



5(b) Therefore there is no firm evidence of asymmetry in 
the MIR emission comparable with that seen in radio and 
X-ray wavelengths. 

4. Discussion 

Very probably the detected MIR flux densities are domi- 
nated by emission from dust grains. 

The contribution of synchrotron radiation can be es- 
timated by extrapolating the radio spectrum given by 
Gaensler et al. l|TM7jl to 10 /im. Taking into account the 
linear increase of radio emission with time, this yields the 
value ^ 5.1 x 10"'^ mJy, two orders of magnitude less than 
the flux densities measured by ISOCAM. 

We can also exclude free-free emission from the debris 
or the thick inner ring as a substantial contributor to the 
observed MIR flux densities. We have estimated the MIR 
free-free continuum from optical recombination line mea- 
surements using published emissivities and recombination 
coefficients for case B of Martin H1988|l . 

For the free-free emission from the ejected material we 
adopted the H^j ffux 2873 days after outburst (Chugai et 
al. . 1199711 assuming a temperature of 10^ K, which is con- 
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offset in arcsec offset in orcsec 

(a) Derived positions of the IR emission (b) Derived sizes of the IR-emitting region 

Fig. 5. Derived absolute coordinates and sizes of the IR-emitting region of the measurements at 12 /im (LWIO) and 
14.3 urn {LW3), overlaid on an HST image of SN 1987A, made in V-band with the WFPC2 on 6. February 1996 
(showing the thick inner ring and emission from the fast expanding ejecta inside). The point source to the east is an 
unrelated star (after convention Star 3, e.g. Walborn et al., 1993). The bright spot on the ring to the south-west is a 
star in the line of sight. 

(a) The coordinates of the infrared emission, derived on the basis of NTT and HST data, are shown as crosses and triangles, 
respectively. The almost circular ellipses, drawn with a solid and a broken line, respectively, show the statistical errors (~ Cf.'lb 
and ~ 0'.'24, respectively at 12 and 14.3 /im) and the total inaccuracy of the derived positions of the IR-emission. The systematic 
error is not well known, because of the edge position of the star, and is simply assumed to be given by 0'.'42, which is the combined 
uncertainty of the pixel coordinates in the inner ISOCAM field after correcting for distortion (Aussel, 1998J. The arrows indicate 
the orientations of the Z-axis of the ISO-satellite. Also shown is the position of the so-called Hot Spot, where at least since 1995 
the emission in optical wavelengths at the inner boundary of the inner ring started to brighten (Lawrence et al.. 120001 . 

(b) The derived sizes of the emitting regions for the two IR measurements with the LWIO and iW5-filter are drawn as ellipses 
centred on the position of SN 1987A. Each of the two axes of the ellipses correspond to the FWHM of the two dimensional 
Gaussian, determined for the IR-source. The orientation and size of the derived size of the IR-emission at 12 ^m agrees very 
well with the projected ellipse of the inner ring in the HST-picture, which is outlined, using the orientation and sizes given by 
Burrows et al. Ill995t . 



sistent with the measured balmer decrement of 3.8. The 
predicted flux density at 10 /im is roughly three orders of 
magnitudes lower 3.2 x 10^'^ mJy). Even in the case of 
a higher temperature of 10'* K the free-free emission from 
the supernova debris in the IR is negligible (~ 0.027 mJy). 

The predicted free-free emission from the thick inner 
ring is brighter, but still insufficient to account for the 
ISOCAM detections. To estimate its contribution to the 
IR the temperature of the whole inner ring was taken to 
be 10** K as was derived from analysis of emission spectra, 
made with the HST at different epochs from a bright part 
of the inner ring (Lundqvist & Sonneborn.,JJ97^. We used 
the flux of the inner ring in R-band, measured with the 
WFPC2 of the HST 3268 days after outburst (Soderberg 
et al.. 11999(1 . Further, we assumed that 19% of the flux is 



due to the He-line and the rest to the neighbouring Nil 
lines at 6548 and 6583 A, which is roughly consistent with 
the line ratios in the line spectrum of a luminous part of 
the ring presented by Panagia et al. H1996fl . For the flux 
density of the free-free emission at 10 fim from the inner 
thick ring we derived an upper limit of ~ 0.053 mJy. 

The contribution of the synchrotron radiation, the 
free-free emission of the ring and the debris to the IR 
is shown in Fig. |S1 which summarizes the spectral energy 
distribution of the SN 1987A at late times. 

No sensitive observations of the MIR cooling lines at 
late epochs are available. However, according to the the- 
oretical predictions for the evolution of the emission from 
various fine structure lines of the supernova debris (Kozma 
& Fransson, il998bjl their contribution should be negligi- 
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ble. We conclude that the measured IR emission is most 
probably thermal emission from dust. This emission could 
either be due to grains in the shocked gas in the cir- 
cumstellar environment or originate in condensates heated 
mainly by the energy of the radioactive decay of the nu- 
cleosynthesis products. 

To derive dust masses knowledge of the optical prop- 
erties of the grains is needed. The grain emissivity can de- 
pend quite strongly on grain composition and in particular 
in the ISOCAM range on wavelength (see e.g. Draine & 
Lee. lI984|) . Nevertheless, for a first interpretation the dust 
spectrum may be approximated by a modified blackbody 
function of the form already used for the colour correction 
in Sect. 13.11 The flux density of N{a) dust grains with ra- 
dius a at a distance D radiating at a certain temperature 
T is then given by 



(7) 



where the emission coefficient is assumed to be Q\{a) = 
a ^ A"'' with constant factors ^ and /3. The total luminosity 
of such a spectrum is given by: 

L = N{a)'i-Ka^ cjT'^ {Q{a,T)) (8) 

^ r(4)c(4) \-c) 

where fc, h, c and a are the Boltzmann constant, the 
Planck constant, the velocity of light and the Stefan- 
Boltzmann constant. {Q{T,a)) is the Planck-averaged 
emission coefficient and r(a;) and Ci^) the gamma- 
and the zeta- functions. Whereas the emission behaviour 
of spherical small grains at long wavelengths can be de- 
scribed with /? = 2 at shorter wavelengths a better ap- 
proximation is given by /3 = 1 . The fit of a spectrum with 
P = 1 including a colour correction to the flux densities 
measured with ISOCAM is shown in Fig. El The corre- 
sponding temperature is 236.0 ± 11.5 K and the luminos- 
ity of the best fit is L = 2.47 x lO^^ W. As seen in the 
figure the approximation gives a good explanation of the 
observed colours (x^ = 0.14). 

The luminosity can be compared with the rate of de- 
position of energy from the decay of radioactive decay 
products. Assuming the ejecta after 11 years to be op- 
tically thin to gamma rays and taking an initial mass of 
10~^ Mq of '*''Ti which is mainly responsible for the heat- 
ing at late epochs then the luminosity emitted by dust 
grains in the MIR is at least ~ 20% of the deposited en- 
ergy rate L « 1.23 x 10^^ W (Fischera, 200^1. 

The dust mass was estimated from the equation 



M, 



dust 



Lp 



3ar4((Q(T,a)) /a) 



(10) 



where p is the grain density. We used Planck-averaged 
emission coefficients at 236 K of silicate and graphite 
grains with a radius a = 0.01 fim adopted from pub- 
lished values from Laor & Draine l)1993|l which are 
(Qsi) /alum] « 0.42 and (Qcra) /a[nm] « 0.080. The 



grain densities are taken to be 2.3 g/cm"^ for graphite and 
3.2 g/cm'^ for silicate grains. The masses of silicate and 
graphite grains are found to be ~ 1.7 x 10~^ Mq and 
~ 6.8 X 10~^ Mq independent of grain size. 

These values indicate a dust mass which is much lower 
than the minimum dust mass of ~ 10~^ Mq derived for 
the newly formed dust in the metal rich part of supernova 
ejecta CWooden. 11997(1 . It is probable that the bulk of the 
condensates are currently emitting in the FIR/submm and 
were too cold to be visible to ISOCAM. A temperature 
of at least ^ 200 K is needed to explain the colours of 
the ISOCAM measurement. By contrast, even as soon as 
~ 1300 days after outburst, the inferred temperature of 
condensates was only ~ 140 K (Bouchet et al. H1991|l . 
Biermann et al., (1992)). The condensate temperature 
should, due to the expansion of the ejecta and weaker heat- 
ing, have been even lower during the ISOCAM observation 
^ 4000 days after outburst. The bulk of the condensates 
should therefore mainly emit in the FIR. Energetically, it 
is feasible that the observed MIR emission could never- 
theless be attributed to a small fraction of condensates 
residing in clumpy regions rich in radioactive nucleosyn- 
thesis products where the local heating would be stronger. 
However, detailed modeling of this case shows that the 
MIR spectral shape would only be reproduced for a small 
range of filling factors, composition and sizes (Fischera, 

mm . 

The fact that the MIR-emission was resolved by 
ISOCAM strongly suggests that the bulk of the detected 
emission does not arise from condensates. The condensates 
are thought to be restricted to the central metal rich core 
of the ejecta (see e.g. Wooden. I1997|l . that is expanding 
with a velocity of 2000 km/s (Kozma & Fransson. ll998a|l . 
At the time of the ISOCAM observations ~ 4000 days af- 
ter outburst the condensates did not cover a region larger 
than - O'.'IS. 

The derived size of 1.501q'24 arcsec and the orientation 
of the emitting region, with an extension consistent with 
the elliptical projection of the thick inner ring fSect. lT^ . 
both point to a predominantly circumstellar origin for the 
dust emission. The MIR extent is furthermore consistent 
with the extensions of the radio (Gacnsler, I1997|l and X- 
ray (Burrows, 2000) emitting regions. The presence of dust 
in the circumstellar environment has also been concluded 
from observations of scattered light correlated with the 
thick inner ring and the associated nebula (Wampler et 
aL. nMTl Crotts et al. lTMSjl . 

Due to the fact that both the radio and the X-ray emis- 
sion are thought to be associated with the blast wave as 
it propagated into an Hll-region of moderate density inte- 
rior of the thick inner ring (Chevalier & Dwarkadas. ll995|l . 
it is likely that the emission in the MIR arises from cir- 
cumstellar dust heated by shocked gas. The observed dust 
temperature of ~ 200 K is consistent with the predictions 
for coUisionally heated grains in equilibrium with a sur- 
rounding hot gas with a density of several 100 cm^'^ (Fig. 7 
in Dwek. [Tl)87|l . Gas density of this order is indeed inferred 
from the observed X-ray emission in the shocked circum- 
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Stellar medium of SN 1987A (Borkowski et aL. ITMTjl . It 
is therefore to be expected that the MIR emission should 
have a detailed morphology similar to the ring-like struc- 
ture as seen in X-rays (Burrows, I2UU(J|I with the centre 
close to the position of the supernova. The derived lumi- 
nosity in MIR is higher than in X-rays (Hasinger et al., 
I1996|l by a factor of ~ 10 as seen also in Fig. This is in 
agreement with values typically observed in galactic SNRs 
(Fig. 6 in Dwek!T988). 

The presence of dust in the circumstellar environment 
is consistent with a red supergiant phase in the time evolu- 
tion of the supernova progenitor (Woosley, 1988). In par- 
ticular, our interpretation of the IR emission is in agree- 
ment with the hypothesis that the Hll-region is wind 
material from this red supergiant phase as suggested by 
Chevalier & Dwarkadas H1995|l . 

5. Summary 

We have presented measurements made with ISOCAM 
at 6.75, 12 and at 14.3 ^m ~ 11 years after outburst of 
SN 1987A. To characterize the IR source we derived the 
absolute coordinates and analysed the extension of the 
emitting region. 

1. The flux densities of the MIR source associated with 
the supernova at 6.75, 12 and at 14.3 fj,m are 0.496 ± 
25.1% mJy, 2.41 ± 4.3% mJy and 2.82 ± 5.5% mJy 
(respectively 3855.2, 3998.7 and 3937.7 days after out- 
burst) (Sect. El) . 

2. The offsets of the IR emission of the remnant 
of SN 1987A at 12 {LWIO) and at 14.3 fim 
(LWS) from the position of the supernova 
(a = 5h35m27.96 8s, S = -69°16'11'.'09 (J2000), 
Reynolds et al., I1995|l . derived relative to an IR- 
emitting star in the ISOCAM fields (see Sect. I3.3jl are 
found to be: 



Aa = -0'.'04±0'.'15,tat. ±0'.'42,y,., 
h0'.'48±0'.'15,tat. ±0'.'42,y,. 



AS = 



and 

Aa 
AS 



-0'.'77 ± 0'.'24stat. ± 0'.'42 
-0'.'09 ± 0'.'24stat. ± 0'.'42 



sys.i 

sys.- 



(LW3) (11) 



(LWIO) (12) 



Due to the edge position of the star the uncertainty of 
the measurement at 12 ^m may be larger. On the basis 
of the absolute coordinates taken in isolation the MIR 
emission could be associated either with the supernova 
or the circumstellar interaction region seen at radio 
and X-ray wavelengths. 

3. The absolute coordinates of the IR-emitting star de- 
rived using observations of the SN 1987A made with 
HST and NTT arc (with a 1 cr uncertainty of - 0'.'09): 

a = 5h35ml8.418s, 5 = -69°16'30'.'65. 

4. From the observation at 12 /im the size of the IR- 
emitting region, assuming a spherical symmetric Gauss 



(FWHM), consistent with the diameter of 1.66±0.03 
arcsec found by Panagia et al. |1991) for the thick in- 
ner ring. Further analysis of these IR data with an el- 
liptical Gauss function also shows a remarkably good 
agreement of the IR-emitting region with the orien- 
tation and the extension of the projected ring even 
though the result is not very significant. The measured 
extension is a clear indication of a predominantly cir- 
cumstellar origin for the IR emission. In particular the 
emission from condensates is not prominent, because 
the condensates would be restricted after 4000 days to 
a region with a diameter of ~ 0'.'18. 
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Appendix A: Theoretical Resolution function 

To estimate the flux, size and coordinates of the IR- 
source we made use of a theoretical point spread func- 
tion (psf) for the ISO-satellite kindly made available to us 
by Dr. Ralph Siebenmorgen. Images of point sources ob- 
served in narrow band LFF-filters of ISOCAM have been 
shown to be (except for LWl) in good agreement with 
this function, which we will refer to as the monochromatic 
psf i Okumur a. II 998|l . For the broad-band filters the psf 
depends however on the source spectrum. Therefore to 
achieve a high accuracy in the determination of the source 
extension we used a theoretical broad-band psf , given by: 



Jf{x) 



J dXXR{\)J{x,\) 
JdXX R{X) Fx ■ 



(A.l) 



Here R{X) is the spectral transmission of an ISOCAM 
filter, taken from the CIA program and Fx the flux density 
(assumed to have a power law Fx oc A"''). J{x, X) is given 

by 



A Fx 

J{x,X) = -^psf{x,X), 



(A.2) 



function for the source, is found to be 1.50 



+0.21 
^0.24 



arcsec 



where A is the effective area of the telescope and psf the 
monochromatic point spread function scaled to unity at 
a; = 0. J{x, A) was derived on a grid of 2048 x 2048 grid 
points sampled at 1/3" x 1/3" using a FFT technique. We 
used values of 60 cm, 17.4 cm and 2 cm for the diameter 



Fischera et al.: Infrared emission towards SN 1987A 11 years after outburst 



Table A.l. FWHM of the measured and the theoretical 
psf 



Filter 


LWIO 


LW3 


FW HM of cuts through the maximum 




3.981 ± 0.084" 


4.980 ±0.012 


FWHM^^^4'Y 


3.951 ±0.252 


4.980 ± 0.063 




3.945 ± 0.026 


4.956 ± 0.058 


relation meas./theor.'' 


1.009 ±0.022 


1.005 ± 0.005 


FW HM of a Gaussian approximation 




4.071 


4.896 




4.035 


4.866 


relation meas./theor. 


1.009 


1.006 



Values taken from Okumura. 119981 
' The given uncertainties correspond to the difference in the 
two values derived in Y and Z direction. 

" Derived from observations of HIC78527 {LWIO) and 
HIC80331 {LW3). 

Due to values from Okumura. 

of the primary mirror, the borehole and the thickness of 
the braces of the tripod, respectively. 

Finally, we convolved the obtained function J-p{x) with 
the used pixel size and the probability distribution of the 
actual pointing position of the satellite during the obser- 
vation {jitter). 

To measure source extensions smaller than the res- 
olution of the instrument it is important to know how 
accurately the theoretical psf represents the actual psf , 
particularly in the kernel. To this end we compared the 
FWHM of the theoretical psf of the LWIO and the LW3 
filters (for spectral index fc = 4) with the actual FWHM 
derived from observations of point sources. The FWHM 
of cuts in Y and Z (in satellite coordinates) through the 
peak of the observed brightness distribution was found 
to be in agreement to within 1% with the FWHM of 
the theoretical psf . A corresponding comparison between 
Gaussian fits to the observed point source and the theo- 
retical psf also agrees within 1% f Table TKA^ . The asym- 
metry of the true psf in the first ring visible in the images 
of point-like sources fOkumura. I1998|l should have only a 
very minor effect. 

An uncertainty of less than 1% is smaller than the 
effect of varying the spectral shape in the calculation in 
the psf . For Gaussian approximations to the theoretical 
psf we found that the dependence of the FWHM on the 
spectral index for the LW3 and LWIO filters is close to 
linear. In a broad range from k w —4 to « 4 the FWHM 
is approximately given by fFischera. I2()()()|l : 

FVKiJM(fc) = 4'.'349 - O'.'OSl fc (iVFlO), , . 

FWHM{k) ^ 5'.'078 - 0'.'048 fc {LW3). ^ ' 

If the theoretical psf of the calibration star (fc = 4) would 
also be taken for the IR emission from the supernova rem- 
nant (which in reality has fc « 1) the FWHM of the psf 
would be underestimated by ~ 6% and ~ 3% for the LWIO 
and LW3 filters, respectively. This can cause substantial 
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systematic uncertainties in the derived source sizes (Fig. 13 
Sect. 1121). 

Appendix B: Method for determination of the 
probability distribution in source size 

To derive the probability distributions shown in Fig. |31 
we measured the variation of x'^i'^o) function of the 
variance ctq of the symmetrical Gauss function assumed 
for the source. We also left flux density {Fq) and coordi- 
nates (a;o) as free variables. Because the source is smaller 
than the resolution of the telescope, we had to convolve 
the source function with the theoretical psf described in 
the previous section. To transform the variation of into 
a probability distribution we approximated the psf itself 
through a two dimensional symmetrical Gaussian function 
with variance ctj . The value of uj was found through a sim- 
ple x^"fit between the psf and the Gaussian. Measuring 
the most probable size ctqq of the source is then equiva- 
lent to measuring the size of the image with the variance 

So = y^fjq^ ± (Tj. In the case of normal distributed errors 
the probability distribution of s with 

^=\A^+^ (B-l) 
is a Gaussian distribution with the variance as'. 

so that the variation of Ax^{s) is a simple quadratic func- 
tion: 

Ax'(s) = x\s) - x'(so) = - so)^ (B.3) 

where x^(*o) — Xmini^)- The measured variation of 
X^(crQ) can therefore be approximated by: 

^'('^q) = ^ ^o)' ± x'(.o). (B.4) 

The most probable value sq, its uncertainty CTs and the 
corresponding value of X^(so) we derived again through a 
X^-fit. Due to uncertainties in the model and the data the 
value So can also be smaller than cj. 

To get the probability distribution over the physical 
sizes of the source we normalized the distribution with: 

„.£.„,,„ = l(l_e,./(fi^)) ,B^5, 

using the error function: 

erf{y) = ^ r dte-'\ (B.6) 

The confidence levels of the maximum sizes of the source 
are simply given by: 

where s is given by Eq. ljB.l|l . 
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Table B.l. Observation log of the calibration stars 



Filter 


LWIO 


LW3 


target 


HIC78527 


HIC80331 


TDT 


189007 


119036 


observer 


CAM CAL 


CAM CAL 


observation date 


15.3.1996 


24.5.1996 


raster 


6x6 


6x6 


pixel size 


3" 


3" 


step size 


2" 


2" 


read out interval 


0.28 s 


0.28 s 


time between pointings 


60 s 


60 s 



B.l. Accuracy of the method 

The accuracy of the used method to derive the extension 
of the emitting region can be estimated by assuming that 
IR source and psf are given by Gauss functions with vari- 
ances (Tq and a J. An uncertainty Aaj in the variance of 
the psf limits a measurement of the extension, neglect- 
ing uncertainties in the data, to sources that are larger 
than (TQ « ^/2Aajaj. The uncertainties in the FWHM 
of the theoretical psf as stated in table lA.ll gives a reso- 
lution hmit of ~ O'.'e {FWHM) for both filters as wiU be 
verified in the next section. 

B.2. Verification of the method using calibration stars 

To verify the method for measuring the angular size of 
the MIR emission associated with SN 1987A we did two 
independent tests. Firstly, we made a control experiment 
by applying the same method to observations of point- 
like source calibration stars (Sect. IB.2.1|I . Secondly, we 
repeated the procedure on the supernova data, but using 
the image of a calibration star as the psf in place of the 
theoretical psf (see Sect. IB.2"2)l . The results both show 
that the model of the theoretical psf allows a measure- 
ment of the size of a Gaussian source with a precision 
which corresponds to an inaccuracy of less than ~ 1% in 
the FWHM. For the two filters we derived a resolution 
limit of about O'.'G at 95% confidence level. 

As measurements of point-like sources we have chosen 
for LW3 and LWIO filters observations of HI78527 and 
of HIC80331, respectively. As in the case of SN 1987A 
the observation provides a sampling of the source with 
an angular resolution of 1" with a pixel size of 3" (Tab. 
IB.1() . The data were calibrated in an almost identical way 
as the observation of SN 1987A with the LW2 filter. The 
only difference was that the deglitching for the calibration 
star data could be done fully automatically using standard 
procedures due to the benign radiation environment en- 
countered for these observations. From the data we first 
derived four images with a sampling of 1" and calculated 
the final one as a weighted average. 



B.2.1. Extension measurements for the calibration 
stars 

The probability distributions in source size of the calibra- 
tion stars were derived from background-subtracted im- 
ages using theoretical psfs for a spectral index of fc = 4 
appropriate for the Rayleigh-Jeans part of the spectrum. 
For the calculation we used a subimage of 11 x 11 pixels 
centred on the source. 

The result is shown in Fig. IB. II As for the supernova 
data we have taken the uncertainties of the pixel values to 
be equal and assigned a value such that the minimum 
corresponds to a confidence level of 90 %. The most proba- 
ble sizes are found to be much smaller than the FWHM of 
the psf. For both filters we found that the source is smaller 
than ^ 0'.'6 with a confidence level of 95%. This is consis- 
tent with the rough estimate of 0'.'6 given in Sect. IB.ll 

B.2.2. Comparison of the size measurements of the 
MIR emission using observed and theoretical psfs 

No point-like calibration source with the same spectrum 
of the source associated with SN 1987A was observed with 
ISOCAM. However, the effect of interchanging the theo- 
retical psf with the observed psf could still be evaluated. 
This was done by comparing calculations for the size prob- 
ability distributions at 12 /im, respectively using as psf 
the observed image of a star and a theoretical psf for 
a spectral index k — A. To achieve for the image of the 
calibration star the same sampling of 1 /3" as for the the- 
oretical psf wc used a cubic interpolation applied to the 
logarithm of the pixel values. It can be seen from Fig. O 
that a most probable size of 1'.'9 is derived both for the 
measured and the theoretical psf. The larger size derived 
by using the measured image of a calibration star instead 
of the theoretical psf appropriate for the supernova is 
caused by the steeper spectrum. 

Appendix C: Position of the reference star 

To determine the position of the IR-emitting star we used 
archival observations with the WFPC2 of the HST, made 
in V-band in February 1996 and three observations (two 
in R one in V-band) made in 1995 of the supernova with 
the SUSI-instrumcnt [SUperb- Seeing Imager) of the NTT 
{New Technology Telescope). The V-band filter was also 
used to determine the correction of the distortion on the 
WFPC2-fields (Holtzman et al.- HM^ . The measurement 
with the NTT was necessary due to the fact that the su- 
pernova and the IR-emitting star were on different CCDs 
of the observation with the WFPC2, which can give an 
additional error in the astrometry. The astrometry on one 
single CCD of the WFPC2 is accurate to O'.'l. The only 
object, whose coordinates are determined with high accu- 
racy is the SN 1987A itself. To remove the offset in the 
astrometry of the WFPC2 data, we have taken the po- 
sition a = 5h35m27.968s, 6 = -69°16'11'.'09 (J2000) from 
Reynolds et al. H1995|l . which is accurate to ~ 70 x 10~^ 
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Fig. B.l. Result of the size measurement for the point-hke sources HIC80331 (iVF3-filter) and HIC78527 {LWIO- 
fiher). The derived probabihty distributions of the size s (see Ea. lB.l|l of the theoretical image B (Eq.[2J) as a function 
of the FWHAI of the Gaussian source is shown in the left panel. The right panel gives the corresponding confidence 
levels for upper limits on source size derived using Eq. IB. 71 



arcsec. We removed the distortion on the WFPC2-field 
using the polynomial correction given by Holtzman et al. 

To derive the astrometry of the observations made with 
the SUSI instrument, we projected all possible star posi- 
tions, seen on a single CCD of the WFPC2, with a tan- 
gential projection on the SUSI field and compared these 
coordinates (xi^yi) with the star positions {xi,yi) on the 
SUSI field with a non linear x^-?ii: 

i 

The positions of the individual stars we derived fitting a 
two dimensional symmetrical Gauss function to the single 
stars in the fields. 

Wc found no indication of a distortion on the SUSI 
field. The final derivations of the star positions in the ob- 
servations with the SUSI instrument and the WFPC2 were 
accurate by cr = 0'.'025 in the case of the PC-field of the 
WFPC2 and by cr = 0'.'049 in the case of the CCD with the 
star (WF3). This is consistent with the accuracy given for 
the single CCDs of the WFPC2 after correcting for distor- 
tion (cr « 0'.'04, Holtzman et aL. ITMSjl . The single CCDs 
indeed show a small offset of up to ~ O'.'IS. On the other 
hand we found no measurable offset for the two CCDs 
with the supernova and the single IR-emitting star. The 
absolute coordinates of the single IR-emitting star in the 
SUSI field are found to be 

a = 5h 35m 18.42s, 5 = -69° 16' 30'.'65. (C.2) 

The coordinates directly derived from the WFPC2 obser- 
vation only differed by ^ 0'.'04 in right ascension and by 
~ O'.'OS in declination. The position of the star relative to 
the SN 1987A should be accurate to 0'.'05. The uncer- 
tainty of the absolute coordinates should be ~ 0'.'09. 
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Fig. 6. Spectral energy distribution of the SN 1987A 
at late times. The measured ISOCAM fluxes are shown 
as open diamonds. They are colour-corrected due to the 
modified black body spectrum F\ oc B\{T)\^^ with a 
temperature of 236.1 K taken as approximation for the 
circumstellar dust emission (solid line). Upper limits in 
broad band at 25, 60, 100, 170 and 200 ^im (short hori- 
zontal lines) are measurements with ISOPHOT made at 
the same epoch as the ISOCAM observations (Fischera, 
2000). The open circles show measured fltrxes at 1.3 mm 
1289 and 1645 days after outburst (Biermann et al. . 11992( 1. 
The La fluxes 3715 and 3734 days after outburst are taken 
from Sonneborn et al. (1998). Upper limits, derived from 
measurements with the SWS-spectrometer aboard ISO 
and shown as filled circles, are taken from Lundqvist et 
al. (1999). Optical fluxes from the thick inner ring (tri- 
angles) and from the debris (squares), derived from ob- 
servations with the WFPC2 of the HST, are taken from 
Soderberg et al. (1999). The synchrotron spectrum is 
an extrapolation from radio measurements, published by 
Gaensler et al. H1997|l . For the Bremsstrahlung in X-ray 
(solid line) we assumed a temperature of 2 x 10^ K and 
scaled the emission measure ^ 2500 days after outburst 
{EM = (1.4 ± 0.4) X 10^^ cm-3, Hasinger et al. ifTMHjl ) 
with EM cx ^2 06 ^j^g ^jj^g Qf i-jjg ISOCAM observa- 
tions. The free- free emission of the thick inner ring and 
the debris are derived from line fluxes, given by Chugai 
et al. (199T), and the flux in R-band, given by Soderberg 
et al. ((1999]). The lines in IR are Fell cooling lines from 
the ejecta at late times derived from theoretical studies 
for IO-'^Mq of -^^Ti (taken from Lundqvist et aL. ITMUji . 
In particular the line at 26 /Ltm, which is thought to be 
the strongest cooling line at late epochs, is of basic in- 
terest as a tracer of the mass of ^^Ti. Using this line the 
3tT upper limit of 51 mJy of the ISOPHOT broad band 
measurement at 25 /zm corresponds to a 3(7 upper limit of 
- 1.5 X 10-* Mq for the mass of ^''Ti CFischera. IMTUjl . 
The emission from condensates is expected to be mainly 
in the FIR (see text). 



